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We propose a phenomenological field theoretical approach to the chemical etching of a disordered- 
solid. The theory is based on a recently proposed dynamical etching model. Through the introduc- 
tion of a set of Langevin equations for the model evolution, we are able to map the problem into 
a field theory related to isotropic percolation. To the best of the authors knowledge, it constitutes 
the first application of field theory to a problem of chemical dynamics. By using this mapping, 
many of the etching process critical properties are seen to be describable in terms of the perco- 
lation renormalization group fixed point. The emerging field theory has the peculiarity of being 
^^self-organized!\ in the sense that without any parameter fine-tuning, the system develops fractal 
properties up to certain scale controlled solely by the volume, V , of the etching solution. In the 
limit V ~* oo the upper cut-off goes to infinity and the system becomes scale invariant. We present 
also a finite size scaling analysis and discuss the relation of this particular etching mechanism with 
Gradient Percolation. Finally, the possibility of considering this mechanism as a new generic path 
to self-organized criticality is analyzed, with the characteristics of being closely related to a real 
physical system and therefore more directly accessible to experiments. 



FAGS numbers; 64.60Ak, 81.65Gf 



I. INTRODUCTION 

Corrosion of solids is an everyday phenomenon of ev- 
ident practical importance The recent development 
of theoretical tools for the study of disordered systems 
and fractals in the context of statistical mechanics (^-||] 
has triggered an outburst of activity in this subject. 

When an etching solution is put in contact with a dis- 
ordered etchable solid, the solution corrodes the weak 
parts of the solid surface, while the hard, stronger, parts 
stay uncorroded. During this process new regions of the 
solid, both hard and weak, are discovered and come into 
contact with the etching solution. As corrosion proceeds 
the etching power of the solution may diminish: indeed, 
if the etchant is consumed in the reaction, etching be- 
comes more and more unlikely until, finally, the solution 
is so impoverished and the solid surface so hardened that 
the corrosion process is arrested. At that moment all 
solid points in contact with the solution are too hard to 
be etched by the weakened etching solution. One of the 
most interesting aspects of this type of phenomenon is 
that the final solid-liquid interface has, in general, a frac- 
tal geometry, at least up to a certain scale |||-^. This 
is precisely the qualitative phenomenology observed in a 
nice experiment on pit corrosion of aluminum thin films 

Recently, a simple dynamical model of etching, captur- 
ing the aforementioned phenomenology, has been pro- 
posed This model has been studied using both 
computational and analytical tools in ||^ , and from these 
studies strong evidence has been provided that the frac- 



tal properties of the solid surface, once the dynamics has 
stopped, are related to isotropic percolation. In princi- 
ple, this is not an obvious result; in fact, at first sight, 
one could think that the interface should be anisotropic 
as there is a preferential direction in which the solution 
advances by etching the solid. 

The purpose of this paper is to provide further theo- 
retical evidence that indeed the critical behavior of the 
model dynamics is related to isotropic percolation. We 
also extend the previous relation to spatial dimensions 
larger than d — 2. To this aim, we shall first review (sec- 
tion 2) two known percolation models that will be useful 
in the forthcoming discussion: (i) dynamical percolation, 
and (ii) gradient percolation (GP). 

Afterwards (section 3), we will define the dynamical 
etching model |^,^ in a circular (spherical) geometry, and 
will derive a phenomenological field theory for it (section 
4). From the analysis of this field theory the parallelism 
with percolation will be set up in a rather clear way, and 
this will provide further theoretical evidence on the con- 
nection between etching and percolation phenomena. 

The approach presented in this paper will allow us to 
study the system self-organization from a field theoret- 
ical point of view, and to verify that, in certain limit, 
the system is self-driven to the neighborhood of a critical 
point without need of any parameter fine-tuning. This 
is a new path to self-organized criticality [0 as will be 
discussed in the last section. 
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II. TWO PERCOLATION MODELS 

In this section we review two different well-known per- 
colation models that will be useful in the discussion of 
the etching processes under consideration. 

A. Dynamical Percolation 

Dynamical percolation is a model proposed for the 
study of the propagation of epidemics in a population, 
and/or for the analysis of forest fires. It is defined as 
follows ||l2|,|l3|. Let us consider a regular square lattice; 
at each site there is a variable that can be in one of three 
possible states (we borrow the language from epidemiol- 
ogy [0): (i) infected sites, (ii) healthy sites susceptible 
to be infected, and (iii) immune sites (non susceptible 
to be re-infected). At time i = a localized seed of in- 
fected sites is located at the center of an otherwise empty 
(healthy) lattice. The dynamics proceeds in discrete time 
steps either by parallel or by sequential updating as fol- 
lows: at each time-step every infected site can infect a 
(healthy) randomly chosen neighbor with probability p 
or, alternatively, heal and become immune to reinfection 
with complementary probability 1—p. Any system state 
with no infected site is an absorbing configuration, i.e., a 
configuration in which the system can get trapped and 
from which it cannot escape |15 Ifl]. It is clear that de- 



pending on the value of p the epidemics generated by 
the initial infection seed will either spread in the lattice 
(for large values of p) or die out (for small values of p) . 
In all cases, the epidemics will leave behind a cluster of 
(healed) immune sites, infinite or finite respectively for 
the two aforementioned cases. Separating the two pre- 
vious phases, there is a critical value of p, < Pc < 1, 
at which the epidemics propagates marginally, leaving 
behind a fractal cluster of immunized sites. It can be 
shown using field theoretical tools (see below) that this 
is a percolation cluster . In this way we have a dy- 

namical model which at criticality reproduces the (static) 
properties of standard percolation. Needless to say, the 
dynamical properties of the dynamical percolation equa- 
tion, do not correspond to any known property of static 
percolation. 

The dynamical percolation model can be cast into the 
following Langevin equation [l^jr^ (or equivalently into 
a field theory MM): 



9ip(x, t) — ^p(x, t) — ap{x, t) I dt'p{x, t') 

Jo 



V2p(x, i) + v/p(x,t)r?(x,t) 



(1) 



where p, (the "mass" in a field theoretical language" ) and 
a > are constants, p(x, t) an activity field describing 
at a coarse grained level the density of infected sites, and 
r]{x,t) a Gaussian white noise. Note the multiplicative 
nature of the noise, because of which the state p{x, t) — 



defines an absorbing state, i.e. dtp{'x.,t) = 0. Note also 
the presence of a non-Markovian term, that constitutes 
the key difference between this equation and the Reggeon 
field theory, characteristic of many other systems with 
absorbing states. This non-Markovian term stems from 
the existence of immunized sites, of which the system 
keeps indelible memory Jl^ , p^ . 

The field theoretical and renormalization group analy- 
sis of Eq.(|l|) can be found in the literature The crit- 
ical dimension is dc = 6, and the exponents, calculated in 
an epsilon expansion, coincide with the well known val- 
ues for percolation calculated using other techniques . 
Apart from the static exponents, also a dynamical expo- 
nent z can be derived from this analysis of dynamical 
percolation [[l3[. 



B. Gradient Percolation 

Gradient percolation |^ (GP) is defined in the fol- 
lowing way: let us consider a bidimensional rectangular 
lattice of lateral sizes L and h respectively, as shown in 

Fig. g. 



p(x) = 
L 
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FIG. 1. Scliematic representation of tlie Gradient Percola- 
tion model. In this case: L = 8 and h = 5. Grey (white) 
rectangles represent occupied (empty) sites. In darker grey 
we indicate the surface of the connected cluster of occupied 
sites. This surface has fractal dimension Df — 7/4 up to the 
characteristic thickness a ~ Vp"^''^^. 

An occupation probability given by p{x) — 1 ~ x/h is 
assigned to sites in column x; this defines a transversal 
constant gradient, Vp = 1/h for the occupation proba- 
bility. Then, at each lattice site (x, y) a random number 
r{x,y) e [0, 1], extracted from an homogeneous distribu- 
tion, and representing the site reluctance to be occupied. 
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is assigned. All sites with r{x,y) < p{x) are declared 
occupied, while the remainings are empty. At the first 
column, a; = 0, all sites are occupied, while there is zero 
occupancy at the last one, x = h (see Fig. |^). After 
identifying all sites as occupied or empty, one detects 
two connected regions (clusters): one (leftmost) with a 
majority of occupied sites possibly surrounding " lakes " 
of empty sites, and another one (rightmost sea) possibly 
surrounding islands. Separating these two regions there 
is an interface (the frontier of the connected cluster of 
occupied sites; it corresponds to the dark sites in figure 
1). The average position of this interface can be shown 
to be at the square lattice site percolation threshold, Pc 
1^,0. In fact, gradient percolation has been used as a 
computational tool to obtain accurate values of percola- 
tion thresholds in different geometries by identifying the 
average position of the interface in sufhciently large lat- 
tices [Q. In the case that we are considering, the fractal 
dimension of the interface, Df = 7/4, can be identified 
as the hull fractal dimension of the critical percolating 
cluster in a two-dimensional lattice |2|]. There is an up- 
per cut-off up to which this fractal behavior is observed; 
it is fixed by the width a which, in its turn, is deter- 
mined by h, and therefore by Vp. It can be shown using 
percolation theory that 

cr - Wp-"" (2) 

where — ^/Df In order to have a well de- 

fined percolation system, with negligible finite size ef- 
fects, the limit L ^ a has to be used. In this way, the 
length /i, determining the value of ct, is the parameter 
that controls the finite size effects; the "thermodynamic 
limit" corresponds to /i — > oo and L — > oo with both 
limits taken in the proper way . One can also estimate 
the variation of p from on the leftmost to the rightmost 
extremes of the wandering interface, Ap: 

Ap^Vp-"'". (3) 

The identity Ap = a * Vp = a/h imposes the follow- 
ing scaling relation among exponents: = 1 — ao-, and 
therefore. 



Let us observe that gradient percolation can also be 
defined in a circular geometry, in which the gradient 
changes with the radial distance to the origin, and the 
cut-off is determined by the width of the roughly circular 
crown in which the interface is inscribed. 

Summarizing, in this section, we have reviewed two 
well-known percolation models. Dynamical percolation is 
a model that, at its critical point, generates dynamically 
a percolation cluster. On the other hand, gradient per- 
colation is a static model, in which an interface appears 
with the same hull fractal dimension of the percolation 
cluster, but with no intrinsic dynamics defined. 



III. DYNAMICAL ETCHING MODEL 

Having introduced the previous two percolation mod- 
els, we go ahead by reviewing the dynamical etching 
model (DEM) at the focus of our study ||,^ . It is defined 
by the following ingredients (see Figure 2). 

(i) The random solid is mimicked by a two-dimensional 
square lattice of finite linear width L and depth Y {Y can 
be arbitrarily large, or even infinite) . Periodic boundary 
conditions in the finite direction are imposed, leading to 
a cylindrical geometry. 

(ii) A random quenched number S [0, 1] (extracted 
from a uniform distribution), assigned to each solid site 
j, represents the site resistance to etching. 

(iii) The etching solution occupies a fixed volume V 
and is initially in contact with the solid through the bot- 
tom boundary, as depicted in Fig. ||, defining a solid- 
solution interface advancing on average in the upward 
direction. 

The solution contains an initial number 7Vet(0) of dis- 
solved etchant molecules. Its concentration at time t is 
C{t) = Net{t)/V. It is assumed that the etching power 
of the solution p(i) is proportional to C{t). Without loss 
of generality the proportionality constant can be fixed 
to unity. Following [||, we assume that etchant parti- 
cles diffuse infinitely fast in the solution (at least much 
faster than the characteristic time scale of etching) and, 
hence, p{t) is taken as spatially homogeneous, i.e., the 
etching power does not depend on the spatial position in 
the solution. 

At each discrete time-step all solid sites located at the 
solid surface and satisfying < p{t) are dissolved (see 
Fig. 1^), i.e., they are removed from the solid, and a par- 
ticle of etchant is consumed for each dissolved site, re- 
ducing in this manner the total etching power. 
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FIG. 2. Schematic representation of the dynamical etch- 
ing model in "cylindrical" geometry before and after the first 
time-step. At this first time-step the (active)sites in contact 
with the solution are i = 1, 2, 3, 4, 5, but only 2,3,5 have a 
resistance lower than the etching power p(0) and then are 
corroded. At the next time step new sites come in contact 
with the solution (the whole second raw if the solution etches 
also in diagonal direction). The etching power diminishes be- 
cause of the consumption of etchant particles. Consequently 
sites 1, 4 stay uncorroded forever. 
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Calling n{t) the number of dissolved solid sites (or 
equivalently the number of consumed etchant particles) 
at time-step t, and N{t) — X]t'=o"'(^') ^^'^ total number 
of etched solid sites up to time t, one can write 

pit + l)^p{t)-^=piO)-^. (5) 

As p{t + 1) < p{t), a site having endured the etching 
attack at time t will also resist at any time t' > t [ pl] |. 
Furthermore, as a consequence of the corrosion process at 
time t, ■m{t) new solid sites, previously in the solid bulk, 
come into contact with the solution at time t + 1. Note 
that they are the sole candidates for corrosion at the next 
time-step. Finally, since the solution has the possibility 
to detach finite solid islands, the global solid surface is 
composed both by the surfaces of the detached islands, 
and by the set of solid sites separating the solution from 
the bulk. This interface is called the corrosion front. A 
more detailed description of the model phenomenology 
can be found in |^ . Here we simply summarize the main 
features of the corrosion front at the arrest time t f . They 
are well represented by GP with Vp ^ L/V: 

(i) the corrosion front shows fractal features with 
1.75 up to a characteristic scale (front thickness) cr; 

(ii) a^(L/^)-i/^/; 

(iii) Pc - p{tf) ^ (i/^^)-"^ with ap ~ {Df - 1)/Df 
(therefore, in the right thermodynamic limit p{tf) Pc). 

Let us introduce here a slight geometrical modification 
of the DEM which makes more clear the connection to dy- 
namical percolation. Instead of considering a cylindrical 
geometry with the etchant solution invading the cylinder 
from the bottom (as in Fig. 2), we consider a flat infinite 
lattice, in which the volume V of the etching solution is 
poured at time i = at an arbitrarily chosen central site 
as schematically shown in Fig. The volume V of the 
etching solution is constant. Observe that with this ge- 
ometry the model has some clear analogies with dynami- 
cal percolation. The main difference is that, in the spher- 
ical DEM, the control parameter (the corroding or infect- 
ing probability) is not a constant but decreases in time as 
the etching process proceeds. As in cylindrical geometry, 
the dynamics can be roughly divided into two regimes : 
a smooth one when p{t) is much larger than pc, and a crit- 
ical one when p{t) approaches Pc- In the smooth regime, 
fluctuations around the average behavior are small while 
in the critical regime fluctuations dominate the dynamics 
1^. Indeed, at early time-steps, the etching power being 
sufficiently larger than pc, it is simple to show 1^ that the 
corrosion front is an approximate expanding circumfer- 
ence centered at the origin, and the number of new solid 
sites coming into contact with the solution at time t sat- 
isfies the approximate relation m{t) ~ 2 tt R{t), where 
R(t) is the maximal radius reached by the corrosion up 
to time time t. As the etchant power is reduced, the 
corrosion front becomes rougher and rougher, until the 
dynamics is finally arrested (see Fig. ||). 




FIG. 3. Pictorial representation of the dynamical etching 
model in "spherical" geometry. 

Since in the smooth regime m{t) 3> 1, we can write 
n(t) ~ p{t)m{t). Hence, within this approximation, it is 
possible to write down the following equation: 



p{t + 1) ~ p{t) - 



27rR{t)p{t) 
V ' 



(6) 



Since in this regime obviously R{t) cx t, we can write 
Eq. ^ in a differential form as follows: 



dp{t) 



dt 



2 TT t 



Pit) 



whose solution is 



p{t) ~p(0)exp (-TTt^/V). 



(7) 



(8) 



From this, the characteristic time of the dynamics is 
seen to be proportional to VV ; i.e. the etching power of 
the solution reaches the value pc in a time tc proportional 
to W. Moreover, as R{t) oc t, at the cross-over between 
the smooth and the critical regime (i.e. when p(t) ~ Pc) 
the solution reaches a distance Rc — VV from the ori- 
gin. By differentiating this expression we conclude that 
the gradient of values of p at which different sites have 
been corroded in the radial direction is proportional to 
R/V. Finally, as R{tc) tx VV, the gradient Vp at tc is 
proportional to 1/VV. In this way, in analogy with the 
cylindrical case, we expect that the geometrical prop- 
erties of the final corrosion front are well represented by 
GP where the gradient of p is dynamically generated. Re- 
placing Vp with 1 / -s/y, the scaling relations studied for 
gradient percolation can be extended to the present case. 
The previous description is valid only for the smooth 
regime, i.e. up to the time at which p{t) w pc- How- 
ever, since the critical regime is shorter than the smooth 
one, we have t/ ^ R{tf) ~ W, where R{tf) is the aver- 
age radius of the final corrosion front, and the previous 
estimations remain valid. In order to check that also 
in the critical regime the radial gradient of the solution 
etching power is given by R/V, it is sufficient to assume 
that during this regime the corrosion front changes from 
a quite smooth geometry to a rougher one, with a final 
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thickness a. Because of the much shorter duration of the 
critical regime one has cr <C R{tf). In this way, during 
the critical regime the solution etches a number of solid 
sites proportional to aR{tf). Therefore, from Eq. (^), the 
variation of the etching power in this regime in average 
is Ap^ aR{tf )/V. 

In conclusion, we have defined a spherical version of the 
DEM, and seen its connection with gradient percolation: 
given the time-diminution of p, the system generates dy- 
namically a spatial gradient of the values ofp at which the 
different sites were etched. Let us finally emphasize that 
if, after the process is arrested, more etchant solution is 
added then the process continues until it is stopped again 
at a value of p around Pc- In this way the disordered solid 
plays the role of a chemical buffer. In the next section 
we present a more theoretical treatment allowing us to 
draw even more precise connections between DEM and 
percolation theory. 



IV. PHENOMENOLOGICAL FIELD THEORY 

In order to construct a field theoretical description for 
the dynamical etching model, a possibility would be to 
write down the master equation defining the dynamics 
and then (using a Poissonian transformation p^j2^ or 
alternatively a Fock-space formalism |^^), derive a gen- 
erating functional Instead of following that strat- 
egy, we prefer here to present a phenomenological set of 
stochastic Langevin equations describing the model at a 
mesoscopic scale. This direct approach, "a la Landau", 
based on the analysis of the main symmetries and con- 
servation laws of the discrete model, has proven very effi- 
cient in describing many other systems related to perco- 
lation, directed percolation and, in general, systems with 
absorbing states [ p^|jl^ . 

Let us consider the following three different local den- 
sities or coarse-grained fields: 

• s(x, describing the local density of material sus- 
ceptible to be etched at any time after t. In the dis- 
crete model there are two types of sites contributing 
to this density: (i) bulk sohd sites and (ii) "fresh" 
interface sites i.e., solid sites freshly arrived to the 
solid-liquid interface (susceptible to be etched at 
the next time-step). 

• q(x, t) is the local density of passivated and in- 
ert material. In the microscopic model this is the 
density of interface sites having already resisted an 
etching trial, i.e., immune or not-susceptible to be 
corroded at any future time-step. 

• c(x, t) is the local density of corroded and replaced 
by solution sites, i.e., the local density of etchant. 

The mean field equations (rate equations) describing the 
evolution of the averaged mean values of these magni- 
tudes are: 



s{t) = -ac{t)s{t) 

g(i) = a(l-p(t))c(t)s(t) 

c(t) = ap{t)c{t)s{t) 



(9) 



where p{t) is the probability to etch an active site at time 
t, and a is a positive constant. In what follows, and with- 
out loss of generality we fix a = 1 . The interpretation of 
the first equation is: in order for the density of suscepti- 
ble sites to change (decrease) in a region, it is necessary 
to have locally both a non-vanishing density of etchants 
and raw solid material susceptible to be etched. This re- 
stricts the dynamics to active regions, i.e., zones in the 
interface separating the etchable-solid and the solution, 
in which non-vanishing local densities of s and of c coex- 
ist. Moreover, the second and the third relations in Eq. ^ 
express the fact that an active site becomes either a c-site, 
with probability p{t) (the corrosion power at time i), or 
alternatively, after healing, a g-sitc with complementary 
probability l—p{t). Note that, as s -I- c -|- g = 0, the to- 
tal number of sites is conserved during the dynamics. It 
is worth stressing that Eq. (^) captures the fact that sites 
resisting an etching attempt remain un-corroded indefi- 
nitely (as occurs in the microscopic model). In fact, the 
number of g-sites grows monotonously until the etching 
process is arrested. 

Observe that we have written so far mean field equa- 
tions in which spatial dependence and fluctuations are 
not taken under consideration. At this point, it is conve- 
nient to introduce an activity field p(x, f) = c(x, t)s(x, t) 
(or p{t) = c{t)s{t) as long as spatial dependences are 
omitted). 

From Eq. it follows immediately that 



(10) 



In order to implement in our theoretical description the 
diminution of the etching power as the corrosion process 
proceeds, we write p{t) (analogously to Eq.(4)) as 



pit) = Po 



m 

V 



(11) 



where now N{t) = J dx[c(x, t) — c(x, 0)] is the number of 
consumed etchant particles up to time t and, as before, V 
is the solution volume. Integrating in time the equation 
for c in Eq. ^, integrating then in space, and substituting 
the result into Eq. |l^, we obtain 



p{t) = p{0) exp 



/ dt' 
Vjo 



dxV(x',i') 



(12) 



Plugging this result into Eq. ^ and re-introducing the 
spatial dependence of the fields, it is a matter of simple 
algebra to obtain 



atp(x, t) = [p(t)s(x, 0) - c(x, 0)]p(x, t) 



p(x, t) 



dt'(l-p{t'))p{^,t') 



(13) 
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In order to go beyond this mean field description, it is 
necessary to include properly spatial coupling and fluc- 
tuations (as a noise). 

• Spacial coupling 

In principle, the spatial coupling can be taken 
into account by introducing additional terms into 
Eq. |l^. However, because of the isotropic nature 
of the local dynamics, terms not invariant under 
space inversion (as odd derivatives of the fields) 
are not allowed. Also, terms like |Vp(x, t)|" can- 
not appear, given the absence of surface tension in 
the microscopic rules Q. Therefore, only terms as 
V^"p(x, t) and higher powers of them are allowed. 
We introduce in Eq. |l^ the only relevant term in 
the renormalization group optics , namely a dif- 
fusive coupling V^p(x, i), that typically appears in 
continuous descriptions of interacting particle sys- 
tems. It will be checked a posteriori that, as a 
matter of fact, the omitted terms are irrelevant at 
criticality. 

• Noise 

In order to introduce properly the noise term, let us 
consider a small region in which there are k "fresh" 
surface sites in contact with the solution of etch- 
ing power p. Since these fresh sites have random 
independent resistances, the average number of dis- 
solved sites will he p ■ k and the fluctuation from 
this average number will be Poissonian, i.e. of order 
\/p ■ k. This implies that in the continuous descrip- 
tion fluctuations of p(x, t) are proportional to its 
square-root. Consequently, a term p(x, t)ri{x, t) 
has to be added to Eq. |l3|, with ri{x,t) being 
a Gaussian white noise with zero mean and no 
spatio-temporal correlations: < r]{x,t)ri{x' ,t') >= 
^(x — x')(5(t — t'). Deviations from gaussianity and 
higher order corrections can be easily argued to be 
irrelevant in the renormalization group optics, and 
therefore are not taken into account. This type of 
noise, with amplitude proportional to the square- 
root of the activity field, is characteristic of systems 
exhibiting a transition from an active to an ab- 
sorbing phase : let us emphasize that wherever 
the activity field is zero, the dynamics is stopped 

HQ- 

Introducing these two new ingredients in Eq.|l3| one ob- 
tains finally 

dtp{x, t) = [p(i)s(x, 0) - c(x, 0)]p(x, t) 



dt'{l-p{t'))p{x,t') 



VV(x, i) + Vp(x,t)r?(x,t) 



(14) 



up to higher order, irrelevant, terms. 

It is worth stressing that even though the microscopic 
model has originally quenched disorder, it has been pos- 
sible to describe it in terms of a stochastic equation with 



annealed noise. This simplification owes to the fact that 
every site is tested for corrosion at most once. If it sur- 
vives, it will stay un-corroded indefinitely, as previously 
explained. In this way, as every random number is used 
at most once it has not to be stored, there are no time 
correlations in the noise and, consequently, a stochastic 
process with no quenched disorder can be used to cast 
the discrete model. 

It is also important to underline that wherever yo(x, t) 
vanishes, all activity, including fluctuations, ceases in 
Eq.p^. In other words, p(pc,t) = defines an absorb- 
ing state ^6 15 1 . This is just another way of saying that 
in the microscopic model, whenever there is no contact 
between the etchable solid and the solution, i.e., when 
they are separated by an interface of passivated (im- 
mune) solid sites, the dynamics is arrested. Continuous 
descriptions of systems with absorbing states are based 
on equations for the activity. In the neighborhood of any 
absorbing state phase transition the activity is small, and 
series expansions on the activity density (as the ones we 
have used to arrive at Eq.|lJ) are justified JT5| , p^ . 

Observe that, apart from the time dependence, of p(i), 
Eq. |l^ is identical to the Langevin equation describing 
dynamical percolation (see Eq. |l|) ||l^,^. As described 
above, dynamical percolation is a rather well known dy- 
namical process generating percolation clusters with a 
characteristic size determined by p. (see Eq. Let us 
analyze the differences between Eq. ^ and Eq!^. Partic- 
ular attention must be paid to the exponential factor in 
the expression for p{t) (see Eq. |2|) absent in dynamical 
percolation; due to it some of the coefhcients in Eq. n4 
are time dependent, while their counterparts in Eq.^ 
are constants. First we discuss whether this extra time- 
dependence may affect the critical properties. Note that 
p{t) does not fluctuate (as verified in simulations in ||]) 
because it is a smooth function of the integral of the field 
over all the past time and the whole space. Therefore, it 
is a deterministically decreasing time-dependent term; or 
in other words it depends on spatio-temporal integrals of 
the activity field and not on the local activity field itself. 
Hence, this term has not critical fluctuations, and does 
not affect the system critical properties. However, it is 
crucial in order to characterize the temporal crossover 
from the active to the absorbing phase. Indeed, as the 
linear term coefficient in Eq. |lj includes a dependence on 
p{t), for early times it is positive, corresponding to the 
fact that for early times the system is in the supercritical 
regime. As the etching mechanism proceeds, the argu- 
ment of the exponential in Eq. grows in modulus, and 
there is a finite time, tc, at which the linear term coeffi- 
cient takes its critical value. Immediately after, the pro- 
cess becomes subcritical i.e. it reaches the absorb- 
ing phase and, the dynamics tends to be stopped with an 
exponentially fast rate (i.e. the number of active sites de- 
creases exponentially). Furthermore, as more and more 
sites are etched, the linear coefficient in Eq. |lj becomes 
smaller and smaller than its critical value, and the expo- 
nential stopping rate is accelerated till a final time </ , at 
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which an absorbing (blocking) configuration is reached. 
Therefore, the main efltect of the time-dependent hnear 
term coefficient is that, by continuously diminishing, it 
drives the system to the neighborhood of the dynamical- 
percolation-field-thcory critical point, but it does not add 
any new relevant operator that could eventually change 
the universality class. At this point, it is a matter of 
simple algebra to verify that all the terms omitted in our 
derivation are indeed irrelevant at the dynamical perco- 
lation, Eq. |l|, renormalization group fixed point. 

Given the previous discussion, one is allowed to say 
that Eqs. |l^ and |lj define a self-organized dynamical per- 
colation Langevin equation: without tuning any parame- 
ter the dynamics is arrested in the neighborhood of the 
percolation critical point, and critical (fractal) properties 
can be measured up to certain length scale determined 
solely by the parameter V (which controls the decrease 
rate of p{t) ). In the limit V ^ oo the upper cut-off goes 
to infinity; i.e., for sufficiently large values of V the point 
at which the dynamics is stopped occurs at values of p 
arbitrarily closed to its critical value in Eq. |l|. 

Observe that the microscopic process is also self- 
organized: initially p{t) is taken to be super-critical, but 
it decreases monotonously till it reaches a critical value, 
and as soon as the sub-critical regime is reached the pro- 
cess is stopped in an exponential way. Therefore, our 
continuous description reproduces the essential features 
of the microscopic model. 

The long range correlations (generating fractal behav- 
ior) in the dynamical etching process are generated in the 
regime in which the linear coefficient takes values around 
its critical value. As a consequence, it is inferred that 
the fractal properties of the final frozen configuration are 
related to the standard dynamical percolation renormal- 
ization group fixed point in any dimension. 

Up to higher order terms, g(x,i) can be written as 
g(x, t) c>c / dt'p{x,t). This variable, the integral over the 
past history of the activity field, represents the statis- 
tics of immunized sites as described by Janssen , and 
therefore, in our problem the statistics of "surviving" 
solid clusters in is also related to percolation proper- 
ties. Observe that regions of the cluster of corroded sites 
far from the final blocking corrosion front have been cor- 
roded with a value of p larger than its critical value and 
therefore are not critical. For the same reason the fi- 
nal corrosion front can be seen as the external perimeter 
of an invading percolation cluster (the etchant solution) 
with p ~ pc- This explains the value of the fractal di- 
mension Df ~ 1.75 found for this final corrosion front, 
which is nothing but the hull exponent of percolation in 
the lattice geometry under consideration p|,p^. 

Finally, we can perform a finite size scaling analysis of 
our equation in order to determine how different magni- 
tudes scale as a function of the only free parameter, V; in 
particular, we can evaluate the distance from the critical 
point at which the process will be finally stopped. The 
linear coefficient in Eq. |lj is 



/i(x, t) = s(x, O)p(O) exp 
-c(x,0). 



i J'dt' J dx'p(x',i') 



(15) 



For the points where the dynamics is arrested, at tj, it is 
clear that s(x, 0) = 1 and c(x, 0) — 0; i.e. at t = they 
belong to the solid bulk. Hence, for the bulk we have 



fi{t) = p{0) exp 



V /o I dx.' p{x.' ,t') independent of 
X. As usual tc is the time at which the process is critical 
p(tc) = /Zc, and tf the time at which the process is actu- 
ally arrested. For large values of V the argument of the 
exponential in Eq. can be expanded in power series, 
and we have, up to the leading order: 

^(ic) - /u(i/) cx — ^ dt' J dx.' p{x ,t'). (16) 

Now we can use the scaling analysis presented in the first 
part of the paper. In the time interval between tc and 
t f the solution erodes a region of the solid with, roughly 
speaking the shape of a circular hole limited by a crown 
of radius Rc and width a. As argued before Rc oc VV. 
Moreover, since the quantity dt' J dx'p{x',t') gives 
the amount of material tested by the etching solution 
between tc and tf, we can write J^^ dt' J dx'p{x',t') ^ 
Rc a. Therefore p,{tc) — l^itf) oc aj^/V. Writing 
a cx (l/-\/y)^"'' one obtains 



fi{tc)~ fi{tf) cx (VV 



(17) 



The fact that V is finite then implies that cr is finite 
and p{tc) — p^itf) 7^ 0. It is only in the large V limit that 
the process is stopped exactly at the dynamical perco- 
lation criticalpoint. As discussed in section 3, one has 
acr = ^/Df IMiI^^! and therefore the distance from the 
final mass to the critical one scales as Ap, that is as the 
excursion of the occupancy probabilities along the sep- 
aration interface of GP. Therefore we have determined 
not only the universality class, but also established how 
finite size effect operate in the DEM. 

In conclusion, all the critical (fractal) properties of the 
microscopic model can be shown to be related to (dy- 
namical) percolation in any space dimension, by using 
the above presented continuous (field theoretical) repre- 
sentation, and finite size corrections can be evaluated. 



V. CONCLUSIONS 

Summing up, we have found that the dynamical etch- 
ing model proposed by Sapoval et al., is a self-organized 
process describable by a the continuous Langevin equa- 
tion similar to that of dynamical percolation. This 
Langevin equation includes a linear term coefficient that 
decreases monotonously as the etching process goes on. 
In this way as soon as it takes it critical value and enters 
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the "absorbing" regime the etching process is stopped. 
Consequently, the fractal (scale invariant) properties of 
the interface in the etching process are shown to be 
related to the (dynamical) percolation renormalization 
group fixed point. This result is valid in any space di- 
mension. In particular, our analysis permits us to con- 
clude that the etching upper critical dimension is dc = 6, 
as in (dynamical) percolation. On the other hand, we 
have also evaluated the role of finite size corrections that 
are essentially different from those of standard dynamical 
percolation. 

An interesting aspect from a theoretical perspective is 
that the field theory (Langevin equation) describing the 
process is self-organized, in the sense that, without any 
parameter fine-tuning, fractal, scale-invariant properties 
are generated. However, it is only in the limit — ^ oo 
{l/V 0), that the upper cut-off for scahng diverges. 
This is in clear analogy with what occurs in other models 
of self-organization, as sandpiles for which critical be- 
havior is observed in the limit of dissipation and driving 
going to zero . 

The mechanism discussed in this paper constitutes a 
new "path" to self-organized criticality in which 

the control parameter decreases monotonously until it 
reaches the neighborhood of the absorbing-state phase 
transition at which the dynamics is arrested in an expo- 
nential way. This same mechanism will be investigated 
in the context of different types of absorbing-state phase 
transitions (as directed percolation |p^ , |l6| ) in a future 
work. Observe that this scenario has the great advantage 
of being related in a clear-cut way to real physical sys- 
tems, making therefore the observation of self-organized 
criticality much more accessible to experiments. 
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